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Using high resolution angle-resolved photoemission data in conjunction with that from neutron and 
other probes, we show that electron-phonon (el-ph) coupling is strong in cuprates superconductors 
and it plays an important role in pairing. In addition to the strong electron correlation, the inclusion 
of phonons provides a theoretical framework explaining many important phenomena that cannot be 
understood by a strongly correlated electronic model alone. Especially it is indispensable to explain 
the difference among materials. The phonons with the wave number around the (0,q x ) and (q x ,0) 
axes create the d-wave pairing while that near (it, it) are pair breaking. Therefore the half-breathing 
mode of the oxygen motions helps d-wave superconductivity. 

PACS numbers: 79.60.Bm, 73.20.Dx, 74.72.-h 



It has been a long-standing question whether a strongly 
correlated electronic model of the CuC>2 plane, such as 
the t-J or Hubbard model alone, can explain the es- 
sential experimental observation of superconductivity in 
cuprates oxides. On the one hand, such a model has 
been remarkably successful in explaining many impor- 
tant physical properties, most notably the property of 
the undoped insulator and the renormalization of chaeae 
dynamics in it, by spin dynamics from t to J scalefl'EI 
and in predicting a spin gapou. On the other hand, 
important questions have been raised. The first is the 
observation that, while the Cu02 plane conductivity is 
essentially the same for various families otc«prates, their 
T c vary by at least an order of magnitudeBLI. The second 
is the observation that the phonoas and lattice effects are 
clearly present in these materialsH, following the original 
assumption that the Jahn- Teller.-/ JT) polarons might be 
important for superconductivity^. There is currently no 
consensus on the above issues. 

Within the context of cuprates superconductors, it is 
difficult to understand some of the material specific prop- 
erties by considering only the electronic degree of free- 
dom. Fig. la summarizes the systematic in the supercon- 
ducting gap size (A), together with the transition tem- 
perature (T c ) shown in Fig. lc. Unlike T c , which can 
be depressed by phase fluctuations in the underdoped 
rcgimcE^rt 2 ], the superconducting gap essentially reflects 
the pairing strength. It can be clearly seen that the pair- 
ing strength for the p-type cuprates is very strong, with 
gap sizes that are at least an order of magnitude larger 
than conventional superconductors. Further, the maxi- 
mum T c of each family is controlled by the maximum su- 
perconducting gap size. This is most dramatically illus- 
trated by the HgBa2Ca n _iCu n 2 („+i) (Hgl223 for n=3) 
compound, which has a much larger gap size as well as 



T c . Fig la also clearly shows a discrepancy between p- 
and n- type materials, as superconductivity in a n-type 
material appears to be very fragile with a much weaker 
pairing strength. The clcctron-holc asymmetry is partic- 
ularly perplexing in the context of spin pairing only, as 
magnetism is very strong and persists to a much wider 
range in the n-type superconductor!!^]. It is natural to 
ask whether one can gain more insight into this issue by 
considering the lattice degree of freedom. 

Considering cuprates superconductors in a larger con- 
text, there are other reasons to investigate the role of 
phonon as exemplified by the following questions. Why 
the lattice i-effect is so important to the properties of 
manganitealil and nickelatealj, but completely unimpor- 
tant to cuprates? Why all the high-temperature su- 
perconductors with T c larger than 30K, MgF^", doped 
C6otIrcJ, BKBOPJj and cuprates, contains light elements 
of boron, carbon and oxygen and have strong bonds? The 
standard theory of electron-phonon interaction does not 
give high-T c and .cL wave pairing, but the issue of electron- 
phonon couplings was not scrutinized to the same ex- 
tend as the strongly correlated theoretical modelsEil. 

This paper addresses these problems by providing ex- 
perimental evidence pointing toward phonon to be also 
an essential player: (i) the phonon strongly influences 
the electronic dynamics, (ii) the strength of the el-ph 
coupling correlates with that of pairing. These finding 
suggests that both electron-electron (el-el) and el-ph in- 
teractions are essential ingredients for superconductivity 
in cuprates. Especially phonon effect is indispensable to 
explain the difference among materials showing different 
superconducting T c . There are two types of the el-ph in- 
teraction, i.e., the diagonal and off-diagonal one. The for- 
mer is the coupling of the atomic (oxygen) displacement 
to the charge at the copper site, while the other is the 
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FIG. 1. In panel (a) we plot the superconducting sap de- 
termined by photoemission spectroscopji-ifull symbol)tJ and 
tunneling spectroscopy (empty symbols )Ej as reported in the 
literature as well as determined from our data. In panel (b) 
the doping dependence of the coupling constant is reported 
for LSCO, Bi2201 and Bi2212 (right side) and NCCO (left 
side) as determined by angle resolved photoemission, along 
the nodal direction (r — Y). In panel (c) (right side), we 
report the doping dependence of the critical temperature for 
four different p-type systems: Hgl2-33, Bi2212, Bi2201 and 
LSCO as reported in the literatureES. On the left side of 
the sa»e panel, we report the critical temperature for the 
NCCCH The shaded area are a guide to the eyes. 



modulation of the transfer integral due to the phonons. 
Deriving the effective t-J modeLirom the three-band d-p 
model according to Zhang-RiceEHl, we obtain both diago- 
nal and off-diagonal interactions whose ratio are different 
between p- and n-type cuprates. The off-diagonal el-ph 
interaction creates the d-wave pairing while the diagonal 
one tends to suppress it. Therefore the weaker d-wave 
superconductivity in n-type cuprates, where the diago- 
nal el-ph interaction is dominant, is consistent with the 
less phonon contribution to the pairing. We also discuss 
the issue of vertex correction for the electron-phonon cou- 
pling. 

In a strong coupling conventional superconductor the 
superconducting transition temperature is 
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where the notations are standardEa. For T r to increase, 



one in principle should increase the phonon frequency uiq. 
This, unfortunately, leads to two problems. The first is 
the Coulomb repulsion /i* = fi/(l + fihi(EF /ivo)), which 
increases rapidly with increasing u>o, since the retarda- 
tion effect is less effective. The second is that higher 
frequency means smaller A, because A = C/Mujq, where 
C is a constant for a given class of materials. In order to 
have a large A, as well as high u>o, one has to increase the 
constant C, which often leads to structural phase tran- 
sitions. These two factors conspire to limit the value of 
T c in conventional superconductors. There are s&vjepil 
band calculations including the el-ph interactionsEZrEi 
The values shown in Fig. 1 for cuprates make it clear 
that the conventional wisdom with phonons will not. be 
able to generate the observed gap size, and thus TcEj. 

It is evident that eq.(l) can not be applied directly to 
cuprates, where the strong el-el interaction plays essential 
role. However we will use this equation below for a guide 
of the discussion, which is not unreasonable because the 
qualitative dependence of the transition temperature on 
the coupling constant A, phonon frequency lu , and the 
Coulomb interaction /i* remains to be correct. 

Although the conventional el-ph interaction alone can 
not explain the high-Tc, recent high-resolution angle- 
resolved photoemission (ARPES) data suggest that 
phonons may nevertheless be an essential player. By ana- 
lyzing photoemission data from P^S^CaCuOs (Bi2212), 
Bi 2 Sr 2 Cu0 6 (Bi2201) and La^Sr-Cu^ (LSCO), in 
addition to earlier data from Bi2212Ej~L3, we concluded 
that the quasiparticles in p-type cuprates are strongly 
coupled to phonons in the frequency rangejJjO-80meVEa, 
in contrast to that of the n-type materials. We see a 
dramatic change in the " quasiparticle" velocity in the 
form of a break in the dispersion near this energy scale, 
as shown in Fig. 2 (a-c) where the dispersion along the 
nodal direction for three families of p-type cuprates are 
plotted. This result is vew-xobust, and the details of the 
fit can be found elsewhereE^O. In related data, a drop in 
the quasip artic le scattering rate below this energy scale 
is observecfiaH, which is consistent with optics data re- 
porting a capid drop of the scattering rate, 1/t, at sim- 
ilar energy^. The observed behavior, a sudden change 
in the dispersion and a drop in the scattering rate, is 
very reminiscent of quasiparticle coupled to a sharp col- 
lective mode. The famous 41 meV neutron mode and 
phonons are the only two sharp collective modes we are 
aware of. Since this effect is seen in all compounds and 
is not associated with the superconducting phase, as the 
phenomena are seen well above T c , we rule out the neu- 
tron mode possibility^, as neutron modes are only seen 
in YBasCusOr-z (YBCO) and Bi2212 and are only as- 
sociated with superconductivity. We note here that the 
overdoped Bi2201 data was collected at 30K, six times 
higher than T c . Alternatively one may try to explain 
the data by the opening of a gap elsewhere on the Fermi 
surface. However, the superconducting and pseudogap 
in LSCO, Bi2201 and Bi2212 are very different, but the 
"kink" energy is very similar, so this explanation does 
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not work. Thus we are left with phonqns as the only 
surviving candidate to explain our datatJ. 

The phonon interpretation receives additional support 
from comparison between photoemission and neutron 
data. Neutron scattering data from LSCO and YBCO 
suggest that the zone boundary phonon couples most 
strongly to doped chargecZl and softens significantly with 
doping. In the case of LSCO, where a direct compari- 
son between neutron and photoemission is possible, the 
phonon mode appears at 70meV. As shown by the thick 
arrow in Fig. 2c, the phonon energy coincides with the 
"kink" energy in the dispersion, providing a very strong 
piece of direct evidence for the phonon being the mode 
responsible for the effect seen in quasiparticle dynamics. 
In a separate analysis we show the striking similarity be- 
tween photoemission data from Bi2212 with that of the 
Be(0001) surface state where the el-ph coupling is known 
to be stronger. In both cases, a peak-dip-hump structure 
is clearly seen in the spectra, as expected, in theories for 
electrons coupled to a collective modeEj'Ea. This analogy 
again provides strong support for the phonon interpreta- 
tion. In summary, we believe that this ensemble of data 
makes a compelling case for a strong coupling between 
quasiparticles and phonons. 

In conventional theory the phonon correction to the 
electronic velocity is given by v(k) = Vo/(l + A), where vq 
is the bare electron velocity (below the phonon frequency 
cl>o), v(k) is the dressed velocity (above the .phonon fre- 
quency) and A is the el-ph coupling constants. In a real 
experiment, one cannot obtain A without knowing the 
bare velocity that can be different from band structure 
value due to el-el interaction. We have extracted a simi- 
lar quantity which we call A' using the velocity obtained 
from the experimental dispersion above the phonon fre- 
quency to approximate the bare velocity. In this case, 
we define the velocity ratio below and above the phonon 
energy as 1/(1+ A'). Because the high energy velocity is 
an overestimate of the bare velocity, A' is an overestimate 
of A. Within a specific model such as Debye model, one 
can see that this overestimate can be 20-30%. We use 
the experimental quantity A' for this discussion as the 
systematics are the same. We find that the velocity ratio 
of p-type materials is close to two near optimal doping 
and increases with underdoping. While the quantitative 
value is model and procedure dependent, the data sug- 
gest that el-ph coupling is strong in these materials. The 
el-ph coupling constant (A') as a function of doping is 
plotted in Fig. lb. We note that the frequency scale of 
the phonon here is an order of magnitude higher than 
that of conventional superconductors. In principle, given 
the high A' value, it can deliver a pairing strength that 
is an order of magnitude higher with respect to the con- 
ventional superconductors, if one can figure out a way to 
deal with the /x* problem. As we will elaborate later, the 
d-wave pairing helps this issue. 

In contrast to p-type materials, the kink effect is much 
weaker and basically not discernible-along the nodal di- 
rection for n-type materials, Fig. 2dc2l, indicating a much 
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FIG. 2. The quasiparticle dispersions vs the rescaled mo- 
mentum (k') are repotted-, for three p-type materials sys- 
tems; Bi22,12 (panel a,f]M, Bi2201 (panel b)E3 and LSCO 
(panel c)H. The arrow indicates the frequency value ob- 
tained by inelastic neutron diffraction data. The dispersions 
are-pompared with the n-type superconductor NCCO (panel 
d)tj along TY. The rescaled momentum, k', is denned as 
(k — kF)/knomeV- The dotted lines are guides to the eye 
and are obtained fitting the high energy part with a linear 
function. 

weaker el-ph coupling. This significantly weakened effect 
is consistent with optical data, where the dronin 1/r is 
found to be much smaller in n-type materiaoS Intu- 
itively this is not surprising, as these high energy phonons 
must involve in-plane oxygen atoms, as suggested by dif- 
ferent experimental. For p-type materials the doped 
carriers predominately go to the oxygen site, while they 
go to the copper site (upper Hubbard band) in n-type 
material. This qualitative difference is expected to give 
a large difference of the el-ph interaction between n- and 
p-type cuprates as discussed below. 

A glance to Fig. la and Fig. lb makes it clear that there 
is a correlation, as a function of doping, between the 
pairing gap size and the el-ph coupling constant. To 
first order, the most striking feature in both Fig. la and 
Fig. lb is that the electron doped and hole doped mate- 
rials are very different. The much weaker el-ph coupling 
strength explains smaller pairing gap in electron-doped 
material. Within a single band t-J or Hubbard model 
with only nearest neighbor hopping integral, one would 
expect particle-hole symmetry. The other clear correla- 
tion is that the pairing strength of p-type material de- 
creases with doping in Fig. la, as does the el-ph coupling 
constant in Fig. lb. We believe that the correlation seen 
in these figures makes a strong case for the el-ph coupling 
being a key to pairing. 

As for the systematic on the gap size in different p-type 
materials, the situation is more complex with the follow- 
ing factors being potentially important. It is known that 
the material dependent copper and apical oxygen dis- 
tance correlates with T c at optimal doping, T™ ax . It 
has been discussed that this is related to the local sta- 
bility of-the Zhang-Rice singlet and its effective transfer 
integralcj. Another electronic structure analysis reduces 
this down to a correlation_between T™ aa: and the range 
of the intralayer hoppingLJ. Hgl201 is found to have 
not only a larger intralayer hopping than LSCO/Bi2201, 
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but also a larger interlayer hopping. The latter is due 
to the on-top stacking of the CuC>2 layers in Hgl201. 
In LSCO/Bi2201 the interlayer hopping is found to be 
substantially smaller due to the body-centered tetrafe 
onal stacking. Difference in hopping may be a reasoned 
for A in Bi2212/Hgl201 (T™ aa: =90K) and LSCO/Bi2201 
(T™ o:E =30-40K) to segregate into two groups. For the 
much higher A max in Hgl223 there is the additional pos- 
sibility of a resonance between the phonon and the su- 
perconducting gap, being the two energies close to each 
other. Optics experiment saw in fact a drop in 1/r in 
Hgl223 system at ,a-|frequency 30% higher than the one 
observed in Bi2212c3, which scales exactly with the pair- 
ing gap (and T c ) difference. The reason for the differ- 
ence in HgI223 is maybe the reduced strain on the Cu02 
sheets by other layers in this material, a factor related to 
the [A* problem as we will discuss later. 

Now we turn to the more detailed description of the 
theoretical consideration. Fig. 3a shows the displace- 
ment pattern of the zone boundary half-breathing mode 
(q=(7r, 0)), which corresponds to the important phonon 
with frequencies 50-80meV. There are several ways to 
consider its interaction with the electronic system. The 
usual approach is to consider the diagonal el-ph interac- 
tion in the particle number, namely the oxygen displace- 
ment is co uple d to the electron density in the i-th 
Cu orbitaleZlea. When one considers the dispersionless 
phonons with the momentum-independent coupling con- 
stant, this in-plane oxygen displacement is pair-breaking 
for d-wave because it gives rise to the repulsive interac- 
tion between the electrons on the nearest neighbor Cu 
orbitals. This is more appropriate when the carriers are 
doped into the Cu orbital as in n-type cuprates, but not 
when they are doped into the oxygen orbital as in p-type 
cuprates. With the strong el-el interaction, the vertex 
correction to the el-ph interaction becomes crucialE3~Eil. 
Especially this diagonal el-ph coupling is irrelevant in 
the undoped case because the charge fluctuation is com- 
pletely suppressed there. In this case the only relevant 
el-ph interaction is in the spin sector, i.e., the modulation 
of the exchange interaction J due to the phonon j-srhich is 
related to the modulation of £y described belowc3. With 
the small doping concentration x, the diagonal el-ph in- 
teraction is still not so effective because the vertex cor- 
rection gives the reduction factor x to the effective cou- 
pling const antcTLil. Further, the absolute value of the 
dielectric constant e(u>) at the energy of the gap ui = 2A 
is much larger than one for the superconducting materi- 
als, and the diagonal el-ph interaction is expected to be 
screened and reduced considerably. 

Therefore we will focus below on the off-diagonal el- 
ph interaction. For the holes, doped predominantly into 
the oxygen orbitals, the direct el-ph coupling originat- 
ing from the shift of the d-level due to the oxygen dis- 
placement is missing. In this case, the el-ph interaction 
should be derived in the framework of the effective t- 
J model where the strong correlation and the exchange 
interaction are taken into account. A realistic model of 



electronic structure is the d-p model, which considers the 
d- and p- orbitals as well as the strong Coulomb interac- 
tions. When we reduce the d-p model into the effective 
single-band t-J modeL the transfer integral £y between 
the Zhang-Rice (ZR)l23 singlet states, at i and j sites, is 
given by the second order process in the hybridization 
tdp ~ 0.8eV between the d- and p-orbitals as 



where Awn-js the energy difference between the p- and 
d-levelsc3'E3. When the negatively charged oxygen ions 
approach (go away from ) the copper, the d-level energy 
increases (decreases). This reflects the fact that the oxy- 
gen displacement modulates the energy level of d-orbital 
as 

5e d = gu(i) (3) 

where u(%) is the linear combination of the oxygen dis- 
placements surrounding the copper i given by 

u(i) =u x (i - x/2) - u x (i + x/2) + u y (i - y/2) - u y (i + y/2), 

(4) 

where u^(i ± /i/2) (/i = x, y) is the displacement along 
the /x-axis of the oxygen between the two coppers at i 
and i ± fi. On the other hand the energy level of the 
oxygen p-orbitals does not change in linear order in it's. 
Therefore the charge transfer energy A<j p = e p — e p in the 
hole picture is given by Ad p = Ao — gu, which leads to 
the off-diagonal el-ph interaction through eq.(2). Ao is 
about 2eV, while g and u are the el-ph coupling constant 
and the displacement, respectively. Then the modulation 
of the transfer integrals Uj is given by 

***i = -S%(«(<)+«(j'))- (5) 

This is the off-diagonal coupling§l. It is noted here that 
Adp(z) and Ad p (J) ( u(i) and u(j) ) appear symmetrically 
in tij . On the other hand, the modualtions of iy due to 
those of tdp's cancel in linear order in it's. In Fig. 3a, we 
take the other view focusing on the displacement of one 
oxygen, and see the modulation of the transfer integrals 
t^. Taking other sites into account the modulation in 
tij, parallel to the displacement (green colored), cancel 
out and only the tij perpendicular to the displacements 
are modulated. 

It is noted here that the transfer integrals £p P 's be- 
tween the oxygen p-orbitals are also relevant to tij in 
the t-J model, and the modulation of t pp 's due to the 
oxygen displacements gives another channel for the off- 
diagonal el-ph interaction. We believe this contribution 
is the main reason why the off-diagonal el-ph interac- 
tion is different between the p- and n-type materials. 
The wave function of ZR-singlet has more weight on the 
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oxygen p-orbitals in the p-type case and is expected to 
have larger el-ph interaction due to this mechanism. It 
should be noted here that the el-ph interaction corre- 
sponding to the modulation of the exchange interaction 
is also present, antLis expected to be dominant in the 
underdoped regionEj. The exchange interaction can be 
represented as J Si ■ Sj = —JxtjXij + constant in terms 
of the bond variable Xij — So- C\ a Gj a . By replacing 

-JxljXij -►-■/(< Xij > Xij + Xij < Xij >) in the mean 
field approximation, this interaction can be reduced to 
the same interaction as eq.(5). An important feature of 
this exchange modulation el-ph interaction is that it is 
free from the reduction due to the vertex correction by 
the factor of x as shown below in eq.(6). It represents the 
fact that it survives in the half-filled case. In Table 1 are 
summarized the el-ph interactions in p-type and n-type 
materials. 





diagonal el-ph 


off-diagonal el-ph 


p-type 


weak and reduced by x 


strong 


n-type 


strongest but is reduce by x 


weak 



Table 1. Classification of electron-phonon interactions. 



Below we focus on the off-diagonal el-ph interaction, 
which is most relevant in the underdoped region due to 
the vertex corrections and can be represented by the stan- 
dard el-ph Hamiltonian as 

#ci- P h = -7= ^2 g(k, q)u q Cl +qa C k<7 . (6) 

where g(k,q) has the characteristic k, g-dependence 
for off-diagonal interactions as g(k,q) oc (sin(g x /2) + 
sin ( q y 1 2 ) ) x ( cos ( k x + q x ) + cos ( k x ) + cos ( k y + q y ) + cos ( k y ) ) . 
These momentum dependence is essential for the pair- 
creating of d-wave superconductivity as shown below. 

The above off-diagonal el-ph interaction, which is tied 
to the change of bond (rather than site) by lattice vi- 
bration, is significantly enhanced by the strong Coulomb 
interaction on the Cu siteEa. This is expressed by the 
vertex correction for the el-ph interaction, which is rep- 
resented by the following formula for the effective el-ph 
coupling constantEH. 

9e B .(k,q) =l(q)9(k,q) = : ^ ' r (7) 

i - m x (q) 

where n x (g) is the generalized susceptibility for bond 

variable Xij — 53 <r CiaQo, an( i the factor 7(g) = 

— * , . is the usual Stoner enhancement factor in the 

i-,m x (q) 

RPA approximation. As described below, this interaction 
creates d-wave pairing. According to the BCS theory, we 
integrate over the phonons to create the effective el-el 
interaction. 
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FIG. 3. In panel a) we report the displacement pattern of 
zone-boundary half-breathing mode of the CUO2 planes. The 
black dots represent the Cu atoms, while the white ones the 
oxygen atoms. The off-diagonal el-ph interaction are shown. 
The displacement of O at the center modulates the energy 
of d-orbitals at the two Cu sites, Cul and Cu2. The effec- 
tive transfer integral t,j perpendicular to the displacements 
are modulatedEj. In panel b) we report the contour plot of 
the momentum dependence of the pairing e(q). Negative e(q) 
indicates pair-creating region, while positive e(q) indicates 
pair-breaking region. The pairing force is zero (green) around 
the nodal region while it decrease gradually, reaching a mini- 
mum (red) approaching the (tt, 0) region. 

H e s = 51 9cs.{k, q)g cS .{k', -q) < u q u^ q > 

k,k' ,q,(T,(T r 

x Ci+ qa C ka Ci_ qa ,C Wa ,. (8) 

Here we have assumed for simplicity the limit of the high 
phonon frequency, which is justified marginally as fol- 
lows. In the present case, uiq = 50 — 80meV is com- 
parable to J = UOmeV which determines the effective 
band-width. We have found that the qualitative conclu- 
sion does not change when we consider the case of more 
realistic bandwidth determined experimentally. 

In the case of s-wave pairing, the Debye cut-off is criti- 
cal and only below it the effective attractive force is effec- 
tive. Therefore the smaller frequency reduce the pairing 
force and enhances the pair-breaking. For the d-wave 
pairing as will be assumed below, on the other hand, 
the momentum-dependence is more important, and the 
sign of the interaction is not relevant. It is because the 
constant part of the effective interaction does not con- 
tribute due to ^ fc < C/cfC-fci >= 0. We have also 
checked that the retardation effect does not change the 
conclusion qualitatively. Our argument is based on the 
fact that the strong correlation and exchange interaction 
leads to d-wave superconductivity. We calculated the 
energy in the BCS approximation for d-wave supercon- 
ductivity, with the order parameter =< Ck^C-ki >= 
(cosfc,, — cos ky)/ 'y/ cos 2 k x + cos 2 k y appropriate for the 
half-filled case. Then we obtain 

2 

< h cS >bcs =-jf^2 \g( k > q)\ 2 \i(q)\ 2 < > A fe A fc+9 

k,q 

= ^2 <u i u -i> e (9)l7(9)| 2 (9) 
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where we have done only the /c-integral in the last line of 
the above equation. It is noted here that < u q U- q >oc 
u>~ 2 and hence the experimental information can be put 
through this factor. Therefore e(q) is weighted more in 
the momentum region where the strong el-ph coupling is 
observed and the dispersion shows softening, i.e., along 
the axis (q x ,0) ancLiO, q y ). In Fig. 3(b) we show the 
contour plot of e(g)E3. It is noted here that the oxygen 
displacement contributes strongest to the d-wave pair- 
ing in exactly the region of q-space, namely along the 
line q = (q x ,0) and q = (0,q y ), where the strongr-Goyr 
pling to the electrons is observed experimentallyEZlE3. 
It should be noted here that u(i) is given by the lin- 
ear combination of the various phonon modes other than 
the half- breathing mode in general, and also even the 
displacements of the apical oxygen play similar roles to 
give off-diagonal el-ph interaction. On the other hand 
the phonons with the momentum around q = (tt, tt) 
are decoupled from the electrons due to the momentum 
dependence of the off-diagonal el-ph coupling constant 
g(k,q) given below eq.(5). In this momentum region, 
the pair-breaking effect originates in the case of diago- 
nal el-ph coupling. In experiments, the anomalous dis- 
persion of this half-breathing mode has been observed 
and the zone boundary region modes show the softening 
and couple to the electrons more strongly .than the zone 
center region modesE3. In another papert3 it is found 
that el-ph coupling is strong from (0.57r,0) to (tt,0). In 
any.xase the finite wave number phonons play important 
roleO. Integration of this range gives pair creating re- 
sults. Especially it is found that |7(g)| 2 strongly depends 
on q, and enhances the coupling around (0, 7r) and (tt, 0) 
points. This corresponds to the tendency towards the 
spin-Pcicrls (dimer) state or stripe formation, and favors 
the d-wave pairing because the pair-creating momentum 
region is enhanced by this vertex correction. 

Now we comment on the n-type cuprates. Although 
our theory can not conclude definitely that el-ph interac- 
tion is weaker in n-type, the smaller kink in Fig. 2 and, 
the smaller scattering rate 1/t(u>) in optical experimented 
strongly suggests it. In particular, the el-ph coupling in 
n-type cuprates is much weaker along the (tt, 7r)-direction 
(although el-ph coupling, or the kink effect is seen along 
other directions), this is consistent with the weaker off- 
diagonal el-ph interaction in the n-type material. Also 
the tunneling data from NCCO has been interpreted by 
the phonons with lower frequencies, which gives lower T c . 
Therefore we can not expect much from the el-ph inter- 
action for the pairing force. It is not clear if this reduced 
pairing force due to el- ph interaction is enough to ex- 
plain the T c in n-type material. It is possible that the 
antiferromagnetic fluctuations cpntribute to the d-wave 
pairing as discussed extensiveljSCj. Therefore certain 
phonon channels can be compatible with d-wave pairing. 
Here we note that our identification of the phonon struc- 
ture in ARPES by the half-breathing q = (tt, 0) phonon 
is mostly motivated by the neutron, experiment showing 
it to soft significantly with dopingj23. The energy consid- 



eration of the photoemission data is also consistent with 
this mode. However, the experimental uncertainty in the 
photoemission data analysis also allows the possibility of 
the lower energy phonons in the range of 40-50 meV. Part 
of this is related to the uncertainty of the superconduct- 
ing gap size which modifies the kink position. Therefore 
the half-breathing mode in our discussion should be re- 
garded as a linear combination of various modes. This is 
particularly true for the in-plane buckling mode (40-50 
meV). According to a band calculatiorCJ, both the half- 
breathing mode and this buckling mode show d-wave pair 
creating tendency near q ~ (0.57T — tt, 0) and pair break- 
ing near q ~ (tt, tt). Raman experiments hav.e-shown very 
strong softening of this phonon through Tjili. We hope 
that our discussion stimulates more theoretical investiga- 
tions on the el-ph interaction in strongly correlated elec- 
tron systems, as the data in Fig. 1 and Fig. 2 strongly 
suggest phonon to be a key player for superconductivity, 
independent of specific theory. It is probably true that 
many phonons need to be considered. 

Fig. 4 proposes a phase diagram that contains essen- 
tial ingredients for the key physics. At very high energy 
and temperature, there is an energy scale (or cross-over 
line) To, which is related with J scale physics caused 
by strong Coulomb interactions. At lower energy scale 
(or temperature) phonons helps pairing which is allowed 
with compatible magnetism to deliver the pairing, yield- 
ing the intermediate pairing energy scale T p h onon (T p h)- 
This two scales energy scheme naturally explains the 
observation that there are two kinds of pseudogaps in 
underdoped cupxates, as strongly indicated by ARPES 
data near (tt, 0)Ex The high enerp; one (hump) connects 
smoothly to that of the insulatorEJ. The low energy_ane 
(leading edge) is related to the superconducting gapl£3. 

A key issue for the classical phonon theory is the prob- 
lem with fi* . Since J has higher energy (Fig. 4), the anti- 
ferromagnetic interaction dictates that the pairing state 
can only be of d-wave symmetry. This suppresses the s- 
wave pairing instabilities of some phonons while cooper- 
ates with the d-wave pairing processes of other phonons. 
Since the d-wave state has a node at the origin, this sig- 
nificantly reduces the fi* problem or even changes sign, 
because there is no amplitude on the same site, and the 
on-site repulsion is not an issue. We note that while /x* in 
eq. (1) is introduced in the context of s-wave supercon- 
ductor as a pair-breaker, the /i* in the discussion here 
should be either positive or negative depending on the 
model without the el-ph interaction. Another critical role 
of the strong Coulomb interaction is that it significantly 
enhances the off-diagonal el-ph coupling near the quasi- 
resonance of the p- and d- levels of cupratesE-i This is un- 
derstood in the vertex correction in eq.(7) due to the near 
instability of the system towards the spin-Peierls tv-ae or- 
dering and/or the bond-centered stripe formatiorO. As 
we have indicated earlier, the high phonon frequency of 
50-80meV and the strong coupling constant can deliver 
a pairing strength that is an order of magnitude higher 
(and thus high-T c ) if the /i* problem can be significantly 
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FIG. 4. A schematic phase diagram for high temperature 
materials is presented. The black line represent the To line, 
which is particle-hole symmetric in the t-J model. The purple 
shaded area represents the range where the phonons become 
important. The superconducting regime (SC) is indicated by 
yellow shaded area and the antiferromagnetic regime (AF) by 
blue shaded area. 

reduced. 

Another factor to differentiate A max among families of 
compounds is the presence/absence of the static distor- 
tion. Aside from theoryoo, that is based on the aver- 
age structural data, this tendency is re-enforced by local 
structural data. As shown by local structural probes, the 
70mcV zone boundary pbsnon couples with dynamic lo- 
cal structural distortions^! The onset temperatures of 
the distortion (T*) are significantly lower in LSCO, cop*- 
pared to Bi2212 and YBCO, likely related to stripesB. 
This coupling with more static structural distortion in 
LSCO is harmful to the pairing and superconductivity 
explaining the lower value of A max , although ujq is com- 
parable. For a similar reason we do not observe any clear 
change of ujo in the Lai.2sNdo.6Sr .i 2 Cu04 (NjdLSCO) 
system with respect to optimally doped LSCOE3. The 
stripes and distortion in NdLSCO are static, further sup- 
pressing pairing. This may also be a factor for A max to 
be the largest in Hgl223 as the Cu02 planes are the flat- 
test there. 

The above interpretation receives support from the 
large change in the onset temperature T* oi-the distor- 
tion (order of 80K) with isotope substitutionEl The op- 
posite change in the T* respect to the one in the T c 
(T* increases with increasing the mass of the isotope), 
is consistent with the idea that more static lattice dis- 
tortions suppress pairing. The local structural distortion 
is likely to be similar in Bi2201 and this may explain 
why the maximum T c (30K) for the single layer is much 
lower than that of the isostructural single layer of Tl and 
Hg compounds. The difference in the dynamics of local 
structural distortion between different families and differ- 
ent dopings is probably related to thejaispiatch between 
the Cu02 plane and the block layerstlTO. As one intro- 



duces more Cu02 planes in a unit cell, their structures 
are less vulnerable to strains from other layers, and this 
may explain why T c increases with number of layers in 
Bi, Tl and Hg based families of cuprates. 

We now briefly discuss the connection between the 
picture which emerges here and those in the litera- 
ture. First there are several authors who considered 
the interplay, between the electron correlation and el-ph 
interactiorliJaOiZj. One important conclusion there is 
that even starting from the momentum independent el-ph 
coupling, the vertex function due to el-el interaction sup- 
presses the large momentum transfer process, and hence 
the transport A* r << Ao. This is expected to resolve the 
apparent contradiction between the ARPES and trans- 
port data. Namely the el-ph interaction of the order of 
Xtr ~ 1 does not appear in the resistivity while A ~ 1 is 
needed to explain the ARPES data. It has been pointed 
out before that the resistivity data is not incompatible 
with |thc presence of el-ph coupling, contrary to general 
belieO. However these works does not take iato account 
the modulation of the exchange interactions which we 
believe is the most important for underdoped region. One 
clue is the ^-dependence of the coupling constant. The 
usual treatment in terms of the 1/N expansion and/or the 
slave boson method leads to the conclusion that the ef- 
fective coupling constant is proportional to x for small x, 
which is not consistent with the observed x-dependcncc 
shown in Fig. 1E3. On the other hand, when one con- 
sider the off-diagonal coupling described in eq. (7), the 
doping will reduce the Stoner enhancement factor of the 
bond-order, and hence the coupling constant is expected 
to be decreased. We have concentrated on the zero- 
temperature superconducting state, and the transport 
properties are beyond the scope of this paper. At finite 
temperature above T c , the spin-charge separation may 
play an important role, and the charge carriers are not 
the simple electrons. If the momenta of the carriers are 
smaller than the Fermi momentum, the large momentum 
phonon scattering does not contribute so much to the re- 
sistivity. However this issue is a subtle one, and we do not 
go further in detail on it. As for the coupling to the spin 
fluctuation as the origin of the kink in the dispersion and 
the peak-hump-dip structure, we believe the small spec- 
tral weight of the 41meV-resonant peak invalidates that 
scenarico. 

Several authors looked for the mechanism of super- 
conductivity in the interplay, bietsgeen the electron corre- 
lation and el-ph interactionoO'Ei Here it should be 
noted that the recent tunneling experiments give evi- 
dence thatpthe-pseudo-gap and superconducting gap may 
be distinctE9il3, while the low energy part of these gaps 
seems to be identicalpJf these picture is correct, which is 
still an open questionllj, this pseudo-gap play an essential 
roles in the physics of underdoped and optimally doped 
cuprates, which can not be captured by the perturbative 
analysis. It is natural to consider that the strongest in- 
teraction, which is the Coulomb interaction and/or the 
spin exchange interaction, determines the global and high 
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energy physics, while the phonons play important roles 
in the lower energy physics. This is the picture we take in 
this paper, assuming for example that the superconduct- 
ing gap is d-wave symmetry, which can not be derived 
by el-ph interaction only. Obviously, diagrams similar to 
Fig. 4 already exist, except that the data in Fig. la and 
Fig. lb strongly suggest that the intermediate pairing 
line is influenced by phonons. 

There are many scenarios for the high energy dynam- 
ics and pseudo-gap formation. One possibility is the an- 
tiferromagnetic-short range order based on the Hubbard 
type modelsE2lEil. The second one is the spin singlet for- 
mation-|bas^d on the RVB pictureEl in terms of the t-J 
modelsutrEJ. In these cases, the role of the phonons has 
been neglected in most of the rlitesature although some 
authors mentioned it implicitlyBEj. For strong coupling 
theories there is no consensus on the lowest energy state 
in these strongly coupled models, with flux phase, d-wave 
pairing, and anti-fcrromagnctically ordered ground states 
being possible candidates. In the intermediate energy 
scales between the pseudo-gap and superconducting gap, 
it is expected that the state can be regarded as the ther- 
mal and /or the quantum mixture of the possible ground 
states such as the staggered flux ancLd-wave supercon- 
ducting state in SU(2) gauge theoryH. The role of the 
el-ph interaction here is to determine the superconduct- 
ing gap size and transition temperature, which differ sig- 
nificantly from material to material. The third scenario 
is the charge ordering and/or the stripe formation dis- 
cussed intensively^. These models are related to what 
we propose, since the stripe approach also includes two 
energy scales and furthermore the half-breathing mode 
can be regarded as the instantaneous creation of the dy- 
namic stripe (Fig. r3ja). In some theories phonons are ex- 
plicitly consideredEj, while in [Others only the electronic 
degrees of freedom are stressed^. For the cases involving 
the Hubbard-Holstein model, the discussion mainly con- 
centrates on the issue of quantum critical points without 
specifying the role of these particular phononsEll For the 
cases where only the electronic aspects of the stripes are 
considered, the two energy scales correspond to the scales 
when the stripes and pairs are formed respectively^. 
In this case the pairing of the carriers is purely elec- 
tronic, driven by quasi- ID behavior. Independently of 
the specifics, it would be interesting to see the role of 
the electronic stripes on the el-ph coupling considered in 
our case. Naturally, charge ordering couples to lattice 
distortions. The fact that the el-ph coupling is strong 
for phonons with q value from (0.57T, 0) to (tt, 0) could be 
interpreted as the system-tendency to couple with stripes 
of 4a to 2a periodicitiesEB. This issue_nLcoupling to 2a 
periodicity has been discussed beforeou. As shown in 
Fig. 3a, the (tt, 0) phonon has a clear ID character with 
2a periodicity. This could be regarded as the instanta- 
neous creation of the dynamic stripes because the charge 
should be accumulated/diluted alternatively along the 
Cu-0 chain perpendicular to the displacement. Alter- 
natively, one can say that the systems tendency to have 



stripes will promote el-ph coupling with desired q, that 
is pair creating for d-wave state. Therefore the stripe 
theory could be regarded as the strong coupling limit of 
the el-ph interaction. 

Before leaving the subject of theory, we also note 
that there are extensive work on phonons only, with- 
out explicit consideration of the strong Coulomb cor- 
relation, notably the polaronic effect and bipolaron 
condensationEJ. Here we think that the d-wave aspects 
are important because of the fi* problem in conventional 
el-ph coupling theory Whatever the specific relationship 
between our finding and the interesting theoretical ideas 
discussed above may be, the data hint strongly that the 
el-ph interaction must be considered explicitly. At the 
same time, we also stress that the el-ph interaction must 
be considered in the context of strong electronic correla- 
tion, as the t-J model has been remarkably successful to 
describe the electronic structure of cuprates down to J 
scalccEl In this sense the strongly correlated t-J model 
does provide a basis for cuprates physics, with phonons 
add to it in delivering d-wave superconductivity with very 
high-T c . 

In summary, we have shown a comparison of photoe- 
mission experiments in conjunction with neutron, and 
other probes, providing direct evidence for strong el-ph 
coupling being important for pairing. The inclusion of el- 
ph interaction explains many experimental observations 
that cannot be understood by strongly interacting mod- 
els of Cu02 planes alone. 
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